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1. INTRODUCTION 

It has been widely recognized that modelling errors can lead 
to sensitivity problems and even divergence of Kalman-Bucy filter. 

In cases where the nominal plant parameters used in the filter design 
are different from the actual plant parameters, the uncompensated 
mismatched steady state Kalman-Bucy filter exhibits bias errors, Athans [1] 
discussed this problem, and presented a brief survey of the various 
schemes to reduce filter sensitivity. In particular Ref. [1] 
introduced the continuous time compensated Kalman filter, a sub- 
optimal state estimator which can be used to eliminate steady state 
bias errors when it is used in conjunction with the mismatched 
steady state (asymptotic) time-invariant Kalman-Bucy filter. The ap¬ 
proach used relies on the utilization of the residual (innovations) 
process of the mismatched filter to estimate, via a Kalman-Bucy 
filter, the state estimation errors and subsequent improvements of 
the state estimate. The compensated Kalman filter augments the mis¬ 
matched steady state Kalman-Bucy filter by the introduction of ad¬ 
ditional dynamics and feedforward integral compensation channels 0 
Satisfactory results of this compensated Kalman filter in a practical 
design have been reported by Michael and Farrar [2], where it was 
applied to the estimation and control of critical gas turbine engine 
variables. 

This note follows the same philosophy and development as [1], 
for the discrete time case. In section 2, we give definitions and 



- 2 - 

and assumptions as well as the definition of the nominal mismatched steady 
state Kalman-Bucy filter. Section 3 analyzes the errors due to mismatching. 
Section 4 contains the main contribution which deals with the development 
of the compensated filter structure and equations. Section 5 contains the 
discussion of the results. 

The elimination of bias errors is accomplished by having accumulators 
(the analog of integrators) acting upon the residuals. Thus if persistent 
bias errors exist due to model mismatching, the accumulators provide the 
necessary corrections so that the bias errors are removed from the estimates 
in steady state. 

The approach selected was dictated by issues of simplicity of design, 
namely the use of constant gains in the estimator realization, and the avoid¬ 
ance of real time parameter estimation, which requires extensive real time 
calculations. All calculations of the constant filter gains can be carried 
out off-line. 
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2. DEFINITIONS AND ASSUMPTIONS 

Only modelling errors in plant parameters will be considered 
throughout this paper. We believe that this is often the case in 
practice. When there are errors in the statistical parameters of the 
underlying random processes, a different approach is required and the 
results are more complicated. 

2.1 Actual Plant Description 

We assume that the actual plant is an n-th order linear time- 
invariant, stochastic dynamical system with state vector x(t)£ R n , 
constant input vector u. £ R m and noisy measurement vector z_(t) £ R r 
described by 


(State Equation) x(t+l) = A x(t) + B u + £(t) (1) 

(Measurement Equation) z_(t) = C_ x_(t) + 9_(t) (2) 

where A, B., C_ are respectively nxn, nxm, rxn constant matrices-. 

The plant noise j^(t) and the measurement noise 0.(t) are white 
Gaussian stationary processes with the following statistics, which 
are assumed to be known to the designer. 


E(?(t))= £ 
E{£(6)S'(t)}= HS st 
E (9_( t)) = 0 

e{q_(s) 0_* ct)} -a6 st 


for all t 
for all s,t 
for all t 


(3) 

(4) 

(5) 

(S) 


for all s,t 



where 6 . denotes the Kronecker delta, and H, the plant noise 
st — 

intensity matrix, is a constant nxn positive semi-definite symmetric 
matrix, while 0, the measurement noise intensity matrix is a 
constant rxr positive definite symmetric matrix. Moreover, £_(t) 
and d(s) are assumed to be uncorrelated for all s, t, i.e. 

E{£(t)9/ (s)> - 0 (7) 

It is further assumed that (with t^ -► - °°) 

E{x(t 0 )} = 0 (8) 


2.2 Model Description 

The actual plant dynamics are defined by the values of the 
three constant matrices A, E[, £ in equations (1) and (2), and 

the noise statistics defined in equations (3)-(8). We shall assume 

that the actual values of A, ]B and C are not known exactly to the 

designer. Rather, nominal values A , B , C are available to the 

—n —n —n 

designer in addition to perfect knowledge of the noise statistics, and 
the value of the constant input vector u. Thus as far as the designer 
is concerned, the model is given by 

x(t+l) = A x (t) + B u + £(t) 

— —n —* —n — — 


( 9 ) 


II 
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For later development, define the following parameter error 
matrices: 


aa 

A 



as 

A - 

A 




—n 

AB 

<111 

B - 

B 




—n 

AC 

A 




c - 

C 

—n 


( 11 ) 

( 12 ) 

(13) 


2.3 Further Assumptions 

We limit the discussions to constant gain filters, as in [1], 
since they represent one of the most practical uses of Kalman- 
Bucy filters from the application viewpoint, and they can readily 
lead into steady state error and stability analysis. The following 
assumptions are necessary for the derivation of the results. 

1. [A,B] and [A ,B ] are controllable pairs. 

— — —n —n 

2. [A,C] and [A ,C ] are observable pairs. 

— — —n —n 

1/2 -1/2 

3. [A,E ] and [A ,- ] are controllable pairs. 

— — —n — 

4. Both A and A are strictly stable matrices, i_e 

— —n 

all of their eigenvalues lie within the unit 

circle. This also implies that (A-I) ^ and (A - I)"" 1 exist. 

— “ —n 

These assumptions are indeed necessary for the rigorous development 
of a unique, stable, steady state Kalman-Bucy filter. 
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2.4 Definition of the Nominal Mismatched Steady-State Kalman 
Buoy Filter (NMSSKBF) 

Let us suppose that the designer constructs the NMSSKBF on the 
basis of the nominal parameter values available to him and the as¬ 
sumed known statistical parameters. Then the state estimate of xitl, 

a yx 

x (t)e R , generated by such a filter, is given by the following 


misma tched filter dynamics. 


x Ct+1 t) = A x (t) + B u 
—n —n-n —n— 


x Ct+1) = x Ct+1 t) + G r Ct+1) , x Ct ) = 0 

—r> —n —n—n —n u — 


r ttl = 2 Ctl - C x Ct t-11 
—n — —n—n 


where G is a constant nxr filter gain matrix given by 


G = 2 C' (0 + C 2 C')” 1 
—n —n —n — —n —n —n 


C14a). 
(14b). 
0.5 1 


(16) 


and 2^ is the nxn constant, symmetric, positive definite solution 
of the algebraic matrix Riccati equation [3] 

A 2 A' - A 2 C'(C 2 C' + 0) _1 C 2 A' + S = 2 (17) 

—n—n—n n—n n —n—n~n — —rr-n—n — —n 

The block diagram of the NMSSKBF is shown in Figure 1. 

The existence of and are guaranteed by our assumptions. 
Also it is well known that the closed-loop filter matrix 

[A - A G C ] = [A - AEC'(0 + C2C ) _1 C J (18) 

—n —n-n—n —n —n—n—n — —n—n—n —n 

is a strictly stable matrix. 




Nominal 

Measurement Residual 



Figure lj Discrete Time Nominal Mismatched Steady State Kalman-Bucy Filter 
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3. ERR05S DUE TO MISMATCHING 

If the NMSSKBF of Figure 1 is used without further compensation , 

severe inaccuracies arise due to the fact that the nominal plant 

matrices A , B , C are used rather than the actual (but unknown) 

—n —n “H 

matrices A, B_, £. These effects are particularly bothersome because 
bias errors in th'e estimates exist. It is instructive to isolate 
these errors, because the structure of the equation suggests that 
compensation techniques can be used. 


3.1 Estimation Errors and their Dynamics 

Define the state estimation error x^(t) induced by the 
NMSSKBF by 


x (t) = x (t) - x (t 1 1 - 1 ) 

“ii — —n 


From equations (1) , (2) , (14) , (15) , and (19) , one can readily 

deduce that the state estimate error x (t) satisfies the stochastic 

-^1 

difference equation 


x (t+1) = A x ( t ) + B u + £(t) - (A x (t) + B u) (20) 
—n - -— —n—n —n— 


Rearranging, one obtains 
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x (t+1) 



+ AB u 


( 21 ) 


3.2 Mean State Estimation Error 

Equation C21) readily allows one to deduce the effects of model 
mismatching upon the estimation errors; bias effects are introduced. 
To compute these, one simply takes expected values in equation (21) • 
We remark that the expectations are not conditional ones, since the 
filter structure has been fixed. Thus 


e{ 2 (t+l)} = [A - A G C ]E{x (t)} 
—n —n —n-rr-n —n 


s{x (t) } 


-A G AC E{x(t) } + /\Eu 

—n—n — — - 

But from equation (1) r 


( 22 ) 



Thus, if some or all of AA f AB , and Ac are non-zero, then one can 
readily conclude from (22) that bias errors exist. Note that the 
effect of the constant input u accentuates these errors, even for 


stable system. From equations C22) and C24) 
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t-t 

E{S (t)} = [A-AGC ] E{5 (t )} 

—n. —n -ii-n -n u 


(25) 


+ Y\ [A-A GC l^CCAA-AG AC] [ X) A S " T B]+Ab}u 
• —n “n —n n — Trn — 


s=t 


T=t, 


Noting that, in view of equations (8) and (14), 

E{x (t_) > = E{x(t n )} - e{£ (t_) } = 0 
—n o — u — —n o — 

equation (25) becomes 

t 


(26) 


E{x (t)> = t. [A -A G C. ] t "* S { [ AA - A G_ AC ] [ £ A S “ T B]+AB>u 

—n —n -II n n — -n—n - . * — — — 

S t — ^ -7 \ 


(27) 


This is in general non-vanishing, and cannot be evaluated since A, B_, C, 
AA, AB, Ac are not known exactly. In particular, as t gets large, the 
mean steady state estimation errors are approached. 


lim E{x(t)} = —[A—I] ^ B u 
tr*° 

lim E{x (t) } = [A -A_G C - I] -1 { [AA-A G AC] [A-I] _1 B-Ab}u 
— n n ~n—n—n — — —rr-n—--— 

t-x" 


(28) 

(29) 


Both of these are non-zero. Therefore, in the mismatched case, there 
exists a non-zero mean steady-state estimation error. Again, due to 
incomplete information, one cannot compute this mean steady-state 

A 

estimate error so as to add it to the NMSSKBF estimate x (t) to arrive 
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at unbiased estimates. 

3.3 Discussion 

The above development indicates that the NMSSKBF should not be 
used without further modification. It necessitates the use of com¬ 
pensation by a suboptimal estimator, because the truly optimal filter 
(that estimates the unknown elements of A, 8_ and C) is an infinite 
dimensional one [3]. The degree of suboptimality has to be related 
to the extra dynamics that are required to improve the performance 
of the nominal mismatched Kalman-Bucy filter. 

The next section presents a compensation technique, which is 
similar to [1]. The central idea is that it augments the NMSSKBF 
by additional filter dynamics and tries to extract further informa¬ 
tion from the residual (innovation) processes. It has the merits 
that all gains can be pre-computed, and it avoids the complexity of 
non-linear estimation (or the extended Kalman filter algorithm) which 
is not always guaranteed to work properly [3]. Moreover, it compen¬ 
sates the biased mean steady state estimate error without sacrificing 
the accuracy of the state estimate, i.e. increase of RMS errors, which is 
very often the case when other techniques are used, e.g., increased 
artifical plant noise covariance matrix I4J. 
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4. THE DEVELOPMENT OF THE DYNAMIC FILTER COMPENSATOR 

First recall the dynamic equations of the state 

estimation error x (t), and the residual process r (t) (equations 
—n n 

(21), (15)). 


X (t+1) = [A -A G C ]x (t) + C(t)-A G0(t) + AA x(t)-A G AC x(t) 
—n —n -n- n -n —n — —n— — — —n-n— — 


+ AB u 


(30) 


r (t) = z(t) - C x Ct t-1) 
—n — —n —n 


(31) 


Using equations (1),(2) and (19), equation (31) can be written as 


r (t) = C x (t) + AC x(t) + 0 (t) (32) 

—n —n —n - 

Note that the residual process r^(t) is linear in the state estimation 

error x (t) . It is reasonable to attempt using equation (32) as a 

—n 

"measurement equation” to obtain an estimate of x^t), the state 
estimation error. However, the existence of the unknown matrices, Aa, 

AB and AC in equations (30) and (32) prohibits us from solving this as 
a linear estimation problem. Thus certain approximations have to be made. 
In what follows, it is shown how one can form such a linear estima¬ 
tion problem, by making approximations with, reasonable physical inter¬ 


pretation. 
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4.1 Philosophy 

Define the time sequences w(t )e R 1 and v(t)e R r as 

follows: 


w(t) = AA sc(t) + AB u 
v(t) = AC x(t) 

Thus equations (30) and (32) reduce to 

x (t+1) = [A -A G C ]x (t) + w(t)—A G v(t) 
—n —n—n—n—n -n — —n—n— 

r (t) = C x (t) + v(t) + 0(t) 

—n —n —n — — 


(33) 

(34) 


+ £(t)-A G0 (t) 

““ (35) 

(36) 


Now it is apparent that equations (35) and (36) form a linear esti¬ 
mation problem, with correlated plant and measurement noise , provided 
that w(t) and v(t) satisfy linear equations. So the next step is to 
develop simple linear equations for w(t) and v(t) . If one is pri¬ 
marily interested in the development of a steady-state constant gain 
filter, one can proceed as follows. 

From equations (28) and (33), we can deduce that 

lim E{w(t) } = [ AB + AA(I-A) ^ B];u = unknown constant (37) 

tr*» 

Also 

w(t+l) - w(t) = AA[(A-I_)x(t) + B u] + AAg( t) 


(38) 
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Equation (37) implies that w(t) must have a nonzero (but unknown) 
mean steady state, while equation (38) indicates that it must contain 
a zero mean driving noise term. Both of these objectives can be 
satisfied if one selects the dynamics of w(.t) to be of the form 

w(t+l) - w(t) = ]£(t) ; w(t Q ) / 0_ (39) 

where ]£(t) is a zero mean, white, Gaussian stationary process, i.e. 


E<X(t) } = 0 

(40) 

E{y(t)Y'(T)} = rfi tT 

(41) 

r = r* >o 

(42) 


Similarly, for v(t), from equations (28) and (34), one can see that 
lim E{v_(t) } ■ - AC (A—I) u = unknown constant (43) 

and 

v(t+l) - v(t) = Ac [ (A-I) x (t) + B u] + AC £(t) (44) 

This leads via to select the dynamics of v(t) by the equation 

v(t+l) - v(t) = X(t) ; v(t Q ) ? 0. (45) 

vdiere X_(t) is a white, Gaussian, stationary noise process with 


statistics 



(46) 


EU(t)} « 0 

E{X(t)X*(T)} = A6 tT (47) 

A = A' >_ £ 

We remark that if AA is known (or can be estimated on the basis 
of physical considerations) one should select 

Jit) = Aa £(t) (48) 

and set 

r = Aa H Aa * (49) 

Similarly, if Ac is known (or can be estimated) one should select 

X.(t) = AC £(t) (50) 

and set 

A = Ac S AC' (51) 

In view of equations (48) and (50), or physical reasoning, 
one should consider y(t) and X^t) as two correlated processes, 
i.e. 


*<!<«*• < T >) - 


(52) 


Indeed if both Aa and Ac are known or can be estimated, one should set 
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Also £(t) should be correlated with x(t) / £(t) . Let 


E{C(t)Y'(T)} = ^ Y 5 tx 
E{£(t)V(T)} = ^ x <S tx 

In particular 



AA' 


if Aa is known (or can be estimated ), and 


(54) 

(55) 


(56) 



(57) 


if Ac is known (or can be estimated). 

The above development allows one, by making certain ap~ 
proximations, to replace equations (30) and (32) by a linear estimation 
problem with plant dynami csl given by 


X (t+l)l 

—n 


"A -A G C I -A G “ 

—n —n—n—n — -n-n 


rt 

_1 


r I o 0 -A G " 
— — — —n—n 

w(t+l) 

= 

£ 1 £ 


w(t) 

+ 

£ 1 £ £ 

_v(t+l) - 


—i 

H| 

ol 

o| 

_1 


- v(t) - 


o| 

o| 

o| 
_1 


£(t+l) 


F 


£(t) 


D 


I<t) 

X(t) 

A(t) 

£(t) 

£(t) 


( 58 ) 


and measurement equation given by 
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r (t) = [C 0 1] 

—n —n — — 


x (t) 


+ 6 (t) 


Thus a Kalman-Bucy filter can be designed to generate estimates of 
j^(t) , w(t) , and v(t) , denoted respectively by xMt) , w(t) , and 

v(t) , based on past measurements of r^T) t <_ T <_ t. 

Recalling that (see equation (19)) 


x(t) = x (t|t-l) + x (t) (60) 

— —n 

it is clear that an improved predicted estimate x{t[t-l), and updated 
estimate x(_t) of x(_t) can he constructed hy 


xCt t-1! = x Ct t'l) + x Ct t-1) 

— 1 —n 1 —n 1 

xCtl = x (_t|t~l) + x t±). 

— —n —n 


Chlal 


C61h) 


4.2 D etails of Constructing the Compensated Filter 

/V 

Using a steady state Kalman-Bucy filter, the estimates x^(t) ' 
w(t) and v(t) are generated via the following equations 

A 

x (t+l|t) = [A -A G C ]x Ct) + wCtl — A G V Ct) - A G 0(0 + HSHL'l‘ 
—n 1 —n —n—n—n —n — -n~n— —n—n-- 

/N 

[r Ct) - C x (t|t-l)-v(t-l)] (62a) 


x (t+1) = x (t+1 t) + L [r (t+1) - C x (t+1 t) - v(t)] (62b) 

—n —n —x —n —n—n — 


w(t+l) = w(t) + L [r (t+1) - C x (t+l|t) - v(t)] 
—w —n —n—n 1 —■ 


v(t+l) = v(t) + L [r (t+1) - C x (t+1 t) - v(t)] 
“V —n —n—n — 
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where 1^, L^, cure respectively nxr, nxr and rxr constant 

matrices (steady state filter gains ), which will be defined below. 
The initial conditions that one might use are 


o 

II 

< 

o 

4-> 

(65) 

A 

w(t Q ) = 0 

(66) 


l« 

ft 

O 

'W 

II 

o 

(67) 


They are chosen because one can view Aa , AB , Ac as random 

matrices with zero means, the nominal values A , B , C used by 

—n —n —n 

the designer represent the a. priori mean values of A, B^ and Q, 
respectively. 

Now referring back to equations C58) and (59)., define for 
notational convenience the matrices 



-Cy - 


a 


2?Y L i£yX ^ 




0 0 0 



( 68 ) 



0 


0 


( 69 ) 
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Then M plays the role of the composite plant white noise Oj>(t)) 
intensity matrix, and N yields the correlation between the composite 
plant white noise $(t) and the measurement white noise 0_(t). 

Let S denote the steady state predicted error covariance matrix 
associated with the estimation problem defined by equations C58), and 
(59). Then S is a symmetric at least semi-positive definite matrix. 
In addition, it is the positive semi-definite matrix solution to the 
algebraic Riccati equation (take e.g. the dual in [5], 


S — (F - D N 0 _1 H] [S_ - £ H' (0 + H S_ H') -1 H S] [F - D N ©“ 1 H] ' 


+ D[M - N Q/V’ ]D' 


(70) 


One can then compute the filter gain matrices L , L and L that appear in 


(62)-(64) by 


L 

-x 

L 
—w 


= S H' (0 + H S H 1 ) 


-1 



(71) 


4.3 Simplification 

The realization of Section 4.2 can be considerably simplified. 
Define the compensated residual process r(t) by 
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r(t) = z(t) - C x(t[t-l) - v(t-l) 

— — — n — 

From eauation (61) we deduce that 

r(t) = 2 (t) - c x (t|t-l) - C x Ct[t-1L - vCt-11 

— — -n—n 1 —n~n 1 ~ 

Then from equation (15) , this becomes 

A 

r (t) = r (t) - Ci (tlt-11 - vtt-11 

— —n —n—n 1 — 

Equations (61), (62), (14) and (74) yield 

A 

x(t+l|t) = X (t+l|t) + X (t+l[t) 

— 1 -n ' —n ' 


—n 


= A x (t|t-l) + B u + A G r (t) 

t"i ■ 1 ■ ■ - n - r*i—n 


+ (A -A G C )x (t) + w(t) -AG v(t) 
—n -n~ n —n —n — —n—n— 

-AG 0(0 + H S H') -1 r(t) 

—n—n-- — 


- A x(t) + B u + w(t) 
— n~- — n— — 


+ [(AG - AGCL - AGL) - AG 0(0 + HSH’) 
—n—n -n—n—n—x n n v —n**n-- 


From, equation (71) , it is easy to see that 

C L + L = HSH 1 (0 + H S H 1 ) -1 
—n—x —v -—- 

After a little algebra, it can be shown that the last term in 
vanishes for all t. That is, equation (75) reduces to 

x(t+l|t) = A x(t) + B u + w(t) 

Thus the compensated state estimate x(t+l) is given by 

A 

x(t+l) = X (t+l|t) + X (t+1) 

— —n 1 —n 

A 

= x (t+1 It) + x (t+1 It) + L r(t+1) 

—n 1 —n. —x— 

= x(t+l|t) + L^rlt+l) 

In addition, equations (63) and (64) become 
w(t+l) = v)(t) + L £(t+l) 


(72) 


(73) 


(74) 


(75) 


a ]r(t) 


equation (75) 


(76a) 


(76b) 


(76e) 


( 77 ) 



- 21 - 


v(t+l) = v(t) + r(t+l) (78) 

Figure 2 shows the realization of equations (76) to (78)• From 

this representation, one can deduce simpler computational procedures for 

the evaluation of the compensated Kalman filter gain matrices L^, L^. 

To be specific, we can reduce the computations of the augmented matrices 

F, D, H, M r N (equations (58), (59), (68), (69)) and the matrix solution 

S_ in the algebraic Riccati equation (70), and the subsequent evaluation 
of the matrices L , L , and L by equation (71). After some algebraic 

—'X —W “V 

manipulations with equations (58), (59), (68)-(71) and (17), one obtains 
the following: 


F 


D N 


0 



r*n 


0 


L £ 

i o 



0 0 1 


D [M-N Q~ 1 N 1 ] D' 


III) 

%Y 

-?Y 

-€y 

r 

Q, . 
-rx 


%x 

A 


( 79 ) 





Figure 2; Discrete Time Compensated Kalman Filter 
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and equations (70), (71) become 

£ “ Hi [S. “ £ H' (0 + H S H ')” 1 H S]F£ + 


(80) 



«-sa’ (0 + H.SH') 


(81) 


4.4 Discussion 

The structure of the discrete-time compensated Kalman filter, 
illustrated in Figure 2, has certain appealing physical aspects. The 
major structural difference between the uncompensated Kalman filter 
of Figure 1 and of the compensated one of Figure 2, hinges upon the 
addition of distinct accumulator loops driven, through appropriate 
gains, by the compensated residual vector £(t+l). Note that an 
accumulator is the discrete time analog of an integrator in the continuous 
time case. In the absence of any accumulators, the residual vector of the 
the uncompensated Kalman filter would exhibit bias errors. In the com¬ 
pensated filter these residual bias errors are accumulated (integrated! 
with appropriate weightings (the gain matrices and in Figure 2 \ 
in two different ways. The accumulator that generates the sequence 
w(t) corrects for bias errors in the predict cycle, to compensate for 
the errors modelled by the matrices Aa and Ab . This makes sense, because 
w(t) was defined by Eq. (33) in terms of Aa and Ab . The accumulator 
that generates the sequence v(t) corrects for bias errors in the 
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residuals due to modelling errors AC in the measurement equation; 
if AC = £, the accumulator channel that generates v(t) would be absent. 

One can then think of w(t) as being a crude estimate of constant 
but unknown disturbances that are additive to the state equations, and 
of v(t) as a crude estimate of constant but unknown bias errors in the 
measurement equation. 

The improved estimate x^(t+l) generated by the compensated Kalman 

filter (Figure 2) is still instantanously influenced by the compensated 

instantaneous residual r_(t+l) through the feedforward gain L^. The 

nominal dynamics A , B , C are still being utilized and bias corrections 

n —n —n 

via w(t) influence the predicted compensated estimate &(t+l|t) and 
hence the residual. 

Obviously the transient performance of the estimation errors of the 

compensated Kalman filter would hinge upon the numerical values of the 

three gain matrices L , L , and L . These in time would not only de- 

—x —w —v 

pend upon the nominal parameters, but upon the way the intensity matrices 
of the white noise sequences y(t) - see eq. (39) - and _^_(t) - see eq. 

(45) - are selected. It is the authors 1 opinion that the suggested 
guidelines for the covariance selection are reasonable, since for most 
practical problems the designer has a reasonable idea of the worst 
possible modelling errors, exhibited in Aa, AB, and AC, from the nominal 
parameters. It is important to stress that the white noise intensity 
matrices and A_ not only depend upon the modelling errors, but also 
upon the intensity matrix E of the original plant noise £(il f shown by 
eqs. £49) and (51). Furthermore, y_(hl and (Jt) are correlated according 
to eqs. (52) to (57). Thus, the guesswork on the part of the designer is 


minimized. 
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5. FURTHER SIMPLIFICATIONS IN PRACTICAL DESIGNS 

At first glance, the compensated Kalman filter might seem some¬ 
what complicated than that of the NMSSKBF alone. We remark that the 
complexity is directly related to the dimensionality of the two 
auxiliary vectors v(t) and w(t) . However, in most cases of practical 
interest, canonical representations can be used to arrive at the 
smallest number of additional accumulators in the realization of the 
compensated Kalman filter. 


5,1 A Single Input Single Output Example 

To illustrate these implications consider the design 
of the compensated Kalman filter in the case of linear, time- 
invariant, stable single input single output plants. Suppose 
that the actual plant is characterized by its transfer function 


G(z) = 


n-1 n-2 

C Z + C -.2 + ...+O..Z + c, 

n _ n-1 2 _ 1 

n n-1 

z + a z +...+a_z + a- 
n 2 -L 


(82) 


Then there are at most 2n unknowns (the aVs an< 3 c^'s) in the P lant 
description. 

Let a." and c. i=l,...,n,denote the nominal values of the a. 
m in i 


and c_^ , respectively. Then the standard controllable representation 


for the actual plant (equation (82)) is 



- 26 - 



u (t) + £(t) 


( 83 ) 


A 


z(t) 



+ 0(t) 


( 84 ) 


and the nominal values of A , B , C are (in the same representation) 

—n —n 


0 


1 


0 


0 


0 


0 


1 . 0 


( 85 ) 


~ a ln “ a 2n .~ an n 


B = 
—n 



( 86 ) 


= [c, c 
In 2n 


c 

nn 


] 


C 

—n 


( 87 ) 
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Thus 


Aa = 


0 

((n-1) rows) 
5a' 


Ab = 0 


AC = 6c* 


where 


<Sa* 

6c' 


[a r’l„' a 2 -a 2n' 


[ °r°ln’ V 0 *’ 


a -a 
n nn 


c -c 
n nn 


] 

1 


Examine the way the vector w(.t)_ is defined (.equation (.331). 

Since AB = 0, w(t) is given by 


w (t) * Aa x (t) 


o 

<Sa'x(t) 


which means that although w(t) is an nxl vector, it really 
only one uncertain parameter, namely the scalar w n (t). 


w n (t) = 6a'x(t) 


( 88 ) 

(89) 

(90) 

(91) 

(92) 

(93) 

contains 


(94) 
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v(t) in 

this case is also a scalar 



v (t) = 5c'x(t) 

(95) 

Following then the procedures described in Section 4, 

one can see 

that the 

scalars w (t) and v(tl should be modelled as 
n 



w^Ct+l) - w n (t) = Y n (t) 

(96) 


v(t+l) - v(t) = A(t) 

(97) 

with 


E{y n (t)} - o 

(98) 


e{A( t)} = 0 

(99) 

and 


E{Y n (t)Y n (T)} - r m s n 

(100) 


E{X(t)X(x)} = 

(101) 


e(y Ct)X(T)} - By x 6 tT 
n f n 

(102) 


£{§.<« Y n (T)> = 

(103) 


E{£(tUCOl - 

(104) 
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Ideally, one selects 


T = 6a' E 6a 
nn- 

A = 6c 1 E 6c 


n 


V 




n 


%x 


5a 1 dc 

H_ Sa^ 

E 6c 


In practice, since <$a and 6c_ are not known one should use their 
worst possible values in equations (105)-(109) to determine the 
covariances. 

Under these conditions, the equations of the compensated 
Kalman filter become (see equations (76)-(78)) 


x(t+l) = A x (t) + 


B u + L r (t+1) + 

—n— x 


w (t) 
n 


w (t+1) = w (t) + L r (t+1) 
n n wn 

v (t+1) = v(t) +L v r(t+1) 


(105) 

(106) 

(107) 

(108) 
(109) 


( 110 ) 

( 111 ) 


( 112 ) 
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with 

r(t) = z(t) - C x(.tjt-ll - v(Jt—1). (113) 

—n— 

where 

a) is a constant nxl gain vector; 

b) L is a constant scalar gain; 

wn 'o • 

c) is a constant scalar gain. 

5.2 Extension to Multi-Input Multi-Output Systems 

The above procedures can be easily extended to the multi-input 
multi-output systems. We outline the step by step procedure that 
should be followed. 

Step 1 : Examine the structure of the nxn matrix Aa for arbitrary 
variations of the actual and nominal plant parameters. Let p , 

0<p<n, denote the number of non-zero rows of Aa . Define an nxp 
matrix P_ according to the following rule. Let j=l,2 ,• ••> P index 
the columns of P. 

(a) . Let i=l, j=l 

(b) . If the i-th row of Aa is nonzero , set the j-th column of P_ 

equal to e^ (the natural basis vector in R n ), set j=j+l, 
and go to (c). if the i-th row of Aa is zero, go to (c) . 

(c) . If i=n stop; otherwise set i=i+l and go to (b). 
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Example ; Suppose Aa has the structure (where x denotes a nonzero 
element) 


AA = 


0 0 0 0 

0 x 0 x 

0 0 0 0 

x x x 0 


Then p=2, and the 4x2 P matrix is given by 


0 0 
1 0 
0 0 
0 1 


Step 2: Examine the structure of the rxn matrix Ac for arbitrary 
variations of the actual and nominal plant parameters. Let q, 
0<q<r, denote the number of non-zero rows of Ac . Define an rxq 
matrix £ according to the following rule. Let k=l,2,... r q index 
the columns of ^ 

(a) . Let i=l, k=l 

(b) . if the i-th row of AC is non-zero, set the k-th 

column of £ equal to e^ (the natural basic vector in 
R r ), set k=k+l, and go to ■ If the i-th row of Ac is 
zero, go to (c). 
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(c). If i=r, stop; otherwise, set i=i+l and go to (b) . 


Exanple: Suppose that Ac has the following structure 


Ac = 


x x 0 0 x 
0 Q 0 0 0 


0 0 0 0 0 


Then q=l and the £ matrix is given by 



Step 3 ; Let (t) e R^, ^ (t) £ , such that 


w^t+l) - w x (t) = (t) 

with 

E [Y-l (t) ] = 0 

E^tt)!-^)] =1^ 


(114) 


(115) 

(116) 


Make the identification 
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w(t) = £ W^(t) 

(117) 



Xj^Ct) = P'Y(t) 

(118) 



£x = £.'£ p 

(119) 

Step 

At 

Let v^t) £ R* 5 , X^t) £ R q such that 


with 


v(t+l) - v(t) = X^t) 

(120) 



E[X L (t) ] = 0 

(121) 



Et^ctjx^x)] = 

(122) 

Make 

the 

identification 




v(t) = £ ^(t) 

(123) 



x (t) = 2 'xct) 

(124) 



a l - 2'A 2 

(125) 

Step 

Sj 

Assume that £(t) , (t) and £^(t) are 

pairwise correlated. 


E[5(t)Y|(T) ] = %Y 1 <S tT (126) 

ElS/tjA^T)] = (12?) 

E[x 1 «=>i. 1 C')i*2 Y x 1 6 tT 


(128) 


Make the identification 


n,_ = P 

-<y 1 -Sy- 


(129) 

(130) 

(131) 


5.3 Reduced Dimension Compensated Kalman Filter; Summary 
of Equations 

The following equations are obtained by substitution of 
equations (114) through (131) into (76) to (78). 

Filter Equations 


X(t+111) = A X(t) + B u + P w. ctl 
— 1 —n _ —n—-1 

(132a) 

x (t+111+1) = x(jt+l|t) + L^r Ct+1) 

C132b) 

w. Ct+1) = w Ct). + L r (t+1) 

—1 —1 —wl—1 

C133a) 

(t+1) = v (t> + L^r^Ct+ll 

a33b) 

r Ct+1). = z Ct+1) - C x(,t+lit) - Q v Ct) 

—1 — n X 

(1341 

Filter Gain Equations 



The filter gain matrices L^Cnxr), (pxr), and L^Cqxrl are given 
by 

= Y_ (0 + I^Y H£ i 



(135) 
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where Y is the unique (symmetric) positive definite solution of the 
(n+p+q)-th order algebraic matrix Riccati equation 


I * £ 2 C - “ - -i ( - + -l- + -2 


(136) 


Fy H^, m 2 are respectively (n+p+q)x(n+p+q), rx (n+p+q) , (n+p+q)x 
(n+p+q) constant matrices given by 


—2 = 


P £ 
I £ 
0 I 


(137) 


—2 



[C 0 
—n — 

2) 



— 




B* 

-Sy 

ii 

-^ A i 




-i _ 


(138) 


(139) 


5*4 Discussion 

The compensated Kalman filter equations presented in Section 5*3 
are the recommended ones for practical design. The additional 
workload associated with the computation of the jP and £ matrices is 
justified in terms of the simplification that results in specific 


problems 



-36- 


From a technical viewpoint the essential difference between 
the filters of Section 4.3 and 5.3 is that in the latter case the 
compensated Kalman filter can be shown to possess all the requisite 
controllability and observability conditions which are necessary 
for the existance of a unique positive definite solution to the al¬ 
gebraic matrix Riccati equation and the guaranteed stability of 
the resultant filter. It also pinpoints the additional number 
Cp+q) of accumulators required to stabilize the uncompensated filter. 

The only "guesswork" required by the designer is the selection 
of the intensity matrices r^(pxp) and (qxq) and the correlation 
covariance matrices (nxp), (nxq) and ft , (pxq). All 


-£y 




U U A 


'1 ’’"l 1 1 

of them scale with the intensity matrix £ of the plant noise C.(t). 

This is evident from equations (49), (51), (53)-(55), (119)., and 

(125)-(128), from which we obtain 



P'Aa 5 Aa'p 

(140) 

II 

rH 

<1 

q*A c £ Ac'q 

(141) 

%Yl 

= £ Aa'p 

(142) 


= £ Ac'Q 

(143) 

ft . 
^ x i 

= P'Aa S AC'Q 

(144) 
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Since £ is assumed known, and P and £ are readily evaluated, the 
selection of F, , A , and , only involves the use 

of judicious estimates for AA and AC by the designer. As 
remarked before, these estimates can be usually found from a worst 
case analysis, i*e., the worst probable deviation of each parameter 
from its assumed nominal value. Of course, there is still an ele¬ 
ment of judicious judgement to be made; at least in this design the 
"guesswork'* can be related to the degree of worst deviation of the 
nominal plant from the actual one. 
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6. CONCLUSIONS 

A method for compensating the mismatched constant gain discrete time 
Kalman filter has been presented. The resultant compensated Kalman filter 
is time invariant, and the gains are all computed off-line, with some 
added complexity to the estimator. The compensation consists of the intro¬ 
duction of feed-forward accumulating (integrating) channels between the 
residual process and the uncompensated mismatched filter. 

The compensated Kalman filter has the property that bias errors in 
the state estimates are eliminated asymptotically. More complex algorithms 
have to be used if unbiased estimates are required for all t. 
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